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ABSTRACT 

We model the small-scale clustering of luminous red galaxies (LRGs; 'Masjedi et al."2006') in the Sloan 
Digital Sky Survey (SDSS). Specifically, we use the halo occupation distribution (HOD) formalism to 
model the projected two-point correlation function of LRGs on scales well within the sizes of their host 
halos (0.016/i"iMpc < r < 0.42/i-^Mpc). We start by varying P{N\M), the probability distribution 
that a dark matter halo of mass M contains N LRGs, and assuming that the radial distribution of 
satellite LRGs within halos traces the NFW dark matter density profile. We find that varying P{N\M) 
alone is not sufficient to match the small-scale data. We next allow the concentration of satellite LRG 
galaxies to differ from that of dark matter and find that this is also not sufficient. Finally, we relax the 
assumption of an NFW profile and allow the inner slope of the density profile to vary. We find that 
this model provides a good fit to the data and the resulting value of the slope is —2.17 ± 0.12. The 
radial density profile of satellite LRGs within halos is thus not compatible with that of the underlying 
dark matter, but rather is closer to an isothermal distribution. 

Subject headings: cosmology: theory — galaxies: elliptical and lenticular, cD — galaxies: fundamental 
parameters — large-scale structure of universe — methods: statistical — surveys 



1. INTRODUCTION 

Galaxy clustering provides important insight into the 
physics of galaxy formation, and can be used to probe 
a broad range of physical phenomena on different cos- 
mological scales. On scales smaller than the virial radii 
of the largest halos (-Rvir ^ 1 Mpc) complex dynamical 
processes are occurring, such as dynamical friction, tidal 
forces, and galaxy mergers. These most certainly have 
an effect on the way galaxies cluster, and will give rise to 
features in the correlation functio n £(r) on small sc ales. 

The Sloan Digital Sky Survey (jYork et al.l 120001 ) has 
produced the largest ever spectroscopic sample of lu- 
minous red galaxies (LRGs). The clustering of these 
galaxies has been considered in great detail on both 
intermediate (iZehavi et al.l l2005f ) and very large scales 
(|Eisenstein et al.ll2005r ). However, it is especially difficult 
to measure the correlation function on small scales due 
to "fiber collision" incompleteness, and deblending is- 
sues. Masjedi et al. (2006, hereinafter M06) were able to 
overcome these observational impasses. They corrected 
for fiber collisions by cross-correlating the spectroscopic 
sample with the imaging sample and testing the results 
against mock data sets. They addressed deblending er- 
rors by introducing artificial galaxies of known magni- 
tudes into the data and studying how they are recov- 
ered by the SDSS software. These adjustments allowed 
them to properly measure ^(r) down to a separation of 
r 15ft.~^kpc. The clustering on these sc ales has also 
been measured for lower luminosity galaxies (|Wang et al.l 
[2006; Li fc Whi te 2009), though these studies did not in- 
clude deblending corrections, which M06 showed can be 
significant. 

In this paper, we model this recently measured 
very small-scale LRG c lustering. On in termediate 
scales (0.3 - iOh-^Mpc), IZehavi et ed] ([2005) measured 
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the two-point cor relation function for 35,000 LRGs. 
iZheng et al.l ()2008l ) modeled this data using the halo oc- 
cupa. t ion distribution (HOD; see, e.g., Pe acock fc SmithI 
'2000': IScoccimarro et al.l [200ll; IBerlind fc WeinbCTgil200a 
iCoorav fc ShethI 120021) framework and found a nice fit. 
However, the extrapolation of their best-fit model to 
smaller separations does not agree with the M06 small- 
scale data: it predicts a correlation function that is too 
low (see M06, Fig. 4). Our motivation is to model these 
innermost data points (0.016 — 0.42^,^-'^ Mpc) to see if 
we can find a model that works. The paper is laid out 
as follows. In §[2l we review the M06 measurement. In 
§ [31 we discuss our method for modeling the small-scale 
correlation function. In § [4l we discuss our results in a 
sequential format: in § 14. li we use four free parameters 
from the probability distribution, P{N\M); in S 14.21 we 
introduce the concentration of satellite LRGs as a new 
free parameter; in § 14. 3|, we allow the inner slope of the 
density profile of satellite LRGs to vary. Finally, in § [5l 
we discuss the implications of our results. 

2. DATA 

M06 measured the small-scale (0.016 - 8/i"^Mpc) 
projected two-point correlation function for a volume- 
limited sample of 24,520 luminous red galaxies in the 
SDSS. The luminosity range of LRGs in the sample 
was —23.2 < Mg < —21.2 and the redshift range was 
0.16 < z < 0.36. Measuring the correlation function on 
such small scales is non-trivial, and requires overcoming 
two main observational hurdles: fiber collisions and de- 
blending. 

The SDSS spectroscopic sample is incomplete due to 
the physical size of the fiber-optic cables used to take 
spectra of targeted galaxies. If two galaxies are closer 
than 55 arcsec on the sky, they cannot both get measured 
spectra. These "fiber collisions" thus result in a mini- 
mum possible pair separation of 55 arcsec. Fiber collision 
incompleteness is reduced by the SDSS tiling method, 
which overlaps the spectroscopic plates to achieve con- 
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tinuous sky coverage. However, this still results in ^ 7% 
of targeted galaxies without measured redshifts. Natu- 
rally, this incompleteness affects galaxy clustering mea- 
surements most severely at very small scales. To get 
around this problem, M06 cross-correlated the spectro- 
scopic LRG sample with the entire sample of LRG tar- 
gets in the SDSS imaging. For every LRG from the 
spectroscopic sample, nearby LRG targets (whether or 
not they have an observed spectrum) were considered 
to be at the same redshift as the spectroscopically ob- 
served LRG. This allowed M06 to assign absolute mag- 
nitudes to the LRG targets and thus decide if they made 
it into the volume-limited sample. M06 then statis- 
tically removed the contribution of galaxies that were 
not at the same redshift, but were considered to be by 
the algorithm, by constructing random samples with the 
same redshift distribution as the spectroscopic sample, 
and cross-correlating them with the LRG imaging sam- 
ple. M06 tested this procedure on mock galaxy catalogs 
and found that it successfully recovers the LRG auto- 
correlation function. 

The SDSS photometry of LRG galaxies is biased in 
cases where pairs of LRGs are separated by tens of kpc 
or less. These galaxies have a region of overlap, and the 
light contained in this region needs to be properly dis- 
tributed between the LRG pair. This process is called 
deblending. Since the LRG sample is defined by lumi- 
nosity cuts (—23.2 < Mg < —21.2), any systematic er- 
rors in deblending will lead to incorrect measurements 
of the correlation function. M06 tested the deblending 
method by introducing artificial, overlapping galaxies of 
assigned magnitudes into the SDSS imaging. They then 
deblcnded the images using t he SDSS imaging pipeline 
PHOTO (iLupton et al.ll2001 fl. and discovered that too 
much light was being systematically allocated to the 
fainter galaxy of the pair. In many cases, this pushed 
the fainter galaxy above the luminosity cut and into the 
LRG sample, whereas it should have stayed out. This 
photometric bias thus led to a boost in the correlation 
function on small scales. M06 corrected their correlation 
function for deblending errors using the results of their 
tests with artificial data. 

After applying these corrections and deprojecting 
Wp(rp), M06 found that the LRG correlation function on 
small scales is a continuation of the £(r) oc power law 
found previously for larger scales (jZehavi et al.l 120051 ) . 
Only their smallest scale data point 10/i~-^kpc) shows 
a downturn, which presumably occurs because at such 
a small scale, it is no longer possible to distinguish two 
merging LRGs from each other. In this paper, we use 
M06's measurements of the projected correlation func- 
tion Wp{rp) from 16 to 420/i^^kpc, which consist of eight 
data points. Restricting ourselves to this range guaran- 
tees that we are safely within the 1-halo regime (i.e., all 
LRG pairs are coming from within the same halos). We 
use the full covariance matrix for this data, which M06 
estimated using jackknife resampling. 

We incorporate the measured number density of galax- 
ies as an additional constraint in our modeling. The 
number density of LRGs with Mg < —21.2 is fig = 
9.73 X 10-5/i3Mpc"^ (|Zheng et al.l [200g) and this pro- 
vides us with a ninth data point. We estimate an error 
for this number density using 50 jackknife samples on the 



sky and obtain cr^g = 1.46 x lO^^h^Mpc ^. 

3. METHOD 
3.1. The Halo Occupation Distribution 

The HOD framework characterizes the bias between 
galaxies (of any class) and mass and is completely de- 
fined by (1) the probability distribution P{N\M) that a 
virialized dark matter halo of mass M will host N galax- 
ies, (2) the relative spatial distribution of the galaxies 
and dark matter within their host halo, and (3) the rela- 
tive velocity distributio n between the gala xies and dark 
matter within the halo (jBerlind fc Weinb erg 2002). We 
may neglect (3) in this study due to the fact that Wp{rp) 
is velocity independent. 

The first moment of P{N\M) is the mean number 
of galaxies as a function of halo mass {N)m, and we 
paramet rize it as a sum of a central and a s atellite com- 
pone nt (jKravtsov et all 120041; [Tinkeii I2007t IZheng et all 
|2008[ ). We assume that there is a minimum halo mass 
cut-off (hereinafter referred to as Mmin) below which a 
halo will always be empty and above which a halo will 
always contain at least a single central galaxy. We next 
assume that the mean number of satellite galaxies is a 
power-law function of mass with a low-mass exponential 
cutoff at Mq. The power-law slope is a and the nor- 
malization is Ml , which represents the mass where halos 
contain, on average, a single satellite galaxy. Specifically, 
the average number of galaxies as a function of mass is 
{N)m = 1 + (A^sat)M, where 

(A^sat)Af = exp[-Afo/(Afhalo - M^in)] X (Mhalo/A^l)" 

The correlation function ^(r) depends on the second 
moment of P{N\M) and so specifying {N)m is not suf- 
ficient. We assume that the actual number of satellite 
galaxies in a halo of mass M follows a Poisson dis- 
tribution of mean {Ngut)- Therefore, the second mo- 
ment is (iVsat (A''sat - l))Af = (A^sat)!/- Our choice of 
a Poisson distrib ution is motivated by th eoretical results 
(jKravtsov et al. 2004; Z heng et al.ll2005D . We note, how- 
ever, that our results will not be very sensitive to the 
specific form of the satellite galaxy distribution because 
of the high mass regime that we are considering in this 
study. LRGs live in halos of mass g reater than the non- 
linear mass M* (jZheng et al.ll200"g ). which means that 
they are on the exponentially declining part of the halo 
mass function. As a result, most of our LRG pairs will 
be central-satellite pairs in halos of mass ^ Mi, rather 
than satellite-satellite pairs in higher mass halos. 

In addition to P{N\M), we must also characterize the 
spatial distribution of galaxies within halos. We nat- 
urally assume that the central galaxy sits at the cen- 
ter of its halo. As for satellite galaxies, we assume at 
first that they trace the density distribution of dark 
matter within their halo, but we eventually relax that 
assumption, as we describe in § El We assume that 
halos have dark matter density profiles that are de- 
scribed by the NFW relation p{v) oc {r/rs)~^{l + r/rs)~'^ 

avarro et al. 1997). The NFW scale radius controls 
where the profile transitions from r^^ in the inner parts 
to in the outer parts, but we parametrize it instead 
through the concentration parameter, which is defined 
as the ratio of a halo's virial radius to its scale radius 
(c = Rvir/rs). Finally, we adopt the concentration - 



3 



mass relation ffiven bv lZheng et"an (|2007D for the modi- 
fication of IBullock et"al] ((2OOII) : c = x (^^)-^, 
where cq = 11, is the non- linear mass at the me- 
dian redshift (M, = 1.2 x 10^'^h-^MQ,z = 0.286) of the 
sample for our choice of cosmology, and /3 — .13. 

3.2. The Galaxy Number Density 

We can calculate the galaxy number density for a given 
HOD by integrating over halo mass and weighting the 
abundance of halos by the mean occupation of galaxies 
{N}m: 



where dn/ dM is the diff e rentia l halo mass function. We 
adopt the iWarren et all (|2006D halo mass function with 
the following cosmological model: fim = 0.25, Vl\ = 
0.75, = 0.04, ho = 0.7, erg = 0.8, n, = 1.0. 

3.3. The Galaxy 2-point Gorrelation Function 

In the halo model, the galaxy two-point correlation 
functio n is given as the sum of the "1-halo" and "2-halo" 
terms ()Zhengjl2004 . 
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The 2-halo term is the contribution of galaxy pairs found 
in separate dark matter halos. Therefore, at scales 
smaller than the virial diameter of the smallest halos con- 
sidered, the 1-halo term compl etel y dominates the corre- 
lation function. IZheng et al.l (f2008t ) found that the min- 
imum halo mass for LRGs is Mmin = 10^^'^Mq, which 
corresponds to virial radii of 0.8 — 0.9/i~^Mpc. There- 
fore, it is sufficient to only consider the 1-halo term when 
modeling ^(r) for our small-scale data. 

The 1-halo term of the two-point correlation function 
can be written as an integral over halo mass, where at 
each mass we add the contribution of central-satellite and 
satellite-satellite pairs: 
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where F^sir) and Fss{r) are the pair separation distri- 
butions of central-satellite and satellite-satellite pairs, 
respectively. The pair distribution for central-satellite 
pairs is essentially the same as the density profile of satel- 
lite galaxies, whereas the satellite-satellite pair distribu- 
tion is equivalent to a co nvolution of the de nsity profile 
with itself. We use the ISheth et"aI1 (|2001f ) calculation 
for the convolution of a truncated NFW profile with it- 
self. As we discussed in § 13.11 central-satellite pairs will 
dominate the LRG correlation function due to the large 
halo masses involved. Therefore, the shape of the small- 
scale LRG correlation function should be very close to 
the shape of the satellite galaxy density profile. 

M06 measured the projected correlation function, so 
we must transform our theoretical ^(r) into w-o( r-o) by 
integrating a.long t he line of sight ((Davis fc Peeblej[l983l : 
iZehavi et al1l20M) : 
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Since we need to integrate ^(r) to large radii in order 
to get w-p{rp), we cannot completely ignore the 2-halo 
term. However, ^(r) is a rapidly declining function of r 
and so u'p(rp) at the small scales we are considering is 
not sensitive to variations in the 2-halo te rm. For this 
reason , we use the best-fit 2-halo term from IZheng et al] 
(|2008f ) instead of calculating it explicitly each time we 
vary our HOD parameters. We tested this by shifting the 
amplitude of the 2-halo term by 20% in each direction 
and we found that w-p[r-p) is not appreciably affected. 

3.4. Probing the Parameter Space 

Now that we have measurements of fig and lyp(rp), as 
well as a mechanism to predict these quantities from a 
given set of HOD parameters (e.g., Mmin, Afo) -^1, a), 
we can probe the parameter space and find the region 
that gives a good fit to the data. We use a Markov 
Chain Monte Carlo (MCMC) method for this purpose. 
Specifically, we adopt the Metropolis-Hastings algorithm, 
which works as follows. The free parameters are given 
initial values and is computed for this starting point 
in parameter space. Steps are then chosen for the pa- 
rameters and is computed for this new location. The 
new location is added to the chain if x^g^ < Xo/d 
n < exp[— (x^eii) ^ Xoid)/^]^ where n is a random number 
between and 1. If these conditions are not met, then 
the old location is repeated in the chain. This process 
then repeats until the chain has converged. We test for 
convergence by starting three chains with different ini- 
tial parameter values and checking whether their final 
parameter distribution functions are in agreement. We 
assume uniform priors for all parameters and do not im- 
pose restrictions on their ranges. 

Once a given chain has converged, we can plot the 
histogram of parameter values in the chain. The most 
likely value for each parameter is given by the mean of 
its distribution, and the associated errors are given by 
the extrema of the middle 68.3% of the distribution. The 
best-fit parameters correspond to the combination of pa- 
rameter values for which x^ is a minimum. For mor e 
details on MCMC techniques, see lDunkley et al.l (|2005D . 



4. RESULTS 



4.1. 



Varying P{N\M) 

As discussed in § [H the IZheng et al] (|2008D fit to 
the intermediate-scale LRG correlation function does not 
match the M06 data points when extrapolated inwards. 
We first investigate whether we can achieve a good fit to 
these new M06 small-scale data points while only vary- 
ing the P{N\M) part of the HOD. We vary the following 
parameters (also defined in 13. ip : 



1. M„ 



the minimum halo mass to contain a central 



galaxy. 

2. Mq - the minimum halo mass to contain satellite 
galaxies. 

3. Ml - the halo mass to contain, on average, one 
satellite galaxy. 

4. a - the slope of the power-law relation between the 
mean number of satellite galaxies and halo mass. 
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Fig. 1. — Model fits to the pro j ected correlation function of LRGs. 
Points show the Masi edi et al.l II2006I) uip(rp) measurements along 
with their jacltknifc errors. The four curves show the best-fit mod- 
els for four sets of free parameters. The PNM model (dotted green 
curve) uses 4 free parameters that describe the probability distribu- 
tion P{N\M) - Mmin, Mo, Ml, and a. The PNMC model (dashed 
blue curve) replaces Mg with /gab which is the concentration of 
the density profile of LRG satellites relative to dark matter. The 
PNMG model (dashed- dotted red curve) replaces this with the in- 
ner slope of the density profile 7. The PNMCG model (solid black 
curve) is the same as the PNMG model with the concentration of 
satellite LRGs /g^i added as a fifth free parameter. The top panel 
shows tOp(rp), and the bottom panel shows the residuals from an 
Tp power law. 



We refer to this model as PNM. Our MCMC rapidly con- 
verges and we find a best-fit model with a reduced of 
6.06 (x^ of 30.3 with 5 degrees of freedom). The dotted 
green curves in Figure [T] show this best fit. It is clear 
that the model provides a poor fit to the data, as it de- 
viates downward from a power law on small scales. We 
therefore find that by varying the P{N\M) free param- 
eters alone we are unable to reproduce the innermost 
M06 data points. We have essentially reproduced the 
disc repancy between th e very small-scale M06 data and 
the iZheng et al.l (|2008[ ) modeling, and thus shown that 
including the small-scale points in the fit does not repair 
the discrepancy. 

4.2. Varying the Concentration of Satellite Galaxies 

Since we cannot reproduce the small-scale LRG clus- 
tering by varying the P(^N\M) distribution, we naturally 
set our sights next on the radial distribution of LRG 
satellites within their dark matter halos. As described in 
§ 13.11 we have assumed that these galaxies trace the dark 
matter halo density distribution, which is in turn de- 
scribed by an NFW density profile. The simplest change 
we can make is to allow galaxies to have a different NFW 
concentration than the halos they occupy. We thus intro- 
duce a new free parameter /gai that relates the satellite 
galaxy concentration Cgai to that of the dark matter halo 
c: 

Cgal = /gal X C. (6) 

In the previous 4 parameter PNM model we found that 
Mq was very poorly constrained. In order to keep the 



same number of free parameters, we fix Mq to M„ii,i, 
setting the exponential cut-off for satellite galaxies to 
occur at M,„iii. Therefore, we now vary the following 4 
free parameters: Mmin, Mi, a, and /gai- We refer to this 
model as PNMC. 

We find a best-fit model with a reduced of 2.52 
(x^ of 12.6 with 5 degrees of freedom). The dashed blue 
curves m Figure [1] show this best fit. The PNMC model 
clearly does better than the PNM model in explaining the 
small-scale LRG clustering; however, it still provides a 
poor fit. Wc find that /gai values of ~ 5—10 are preferred, 
showing that LRGs are more concentrated than the dark 
matter for this PNMC model. This makes sense because 
increasing Cgai means that we are adding more satellite 
galaxies towards the center of halos. This forces the scale 
radius inwards and boosts the amplitude of the inner part 
of ^(r). However, while simply moving more galaxies 
towards the center may aid in fitting the very inner most 
2-3 data points, this can result in a poorer fit to the 
outermost data points. In other words, varying Cgai can 
shift Wp(rp) in the rp direction, but it cannot alter its 
shape, which is fundamentally not a power law in the 
case of an NFW satellite profile (see dashed blue curve 
in bottom panel of Fig. [T|). 

4.3. Varying the Density Profile 

Adopting an NFW form for the density profile of satel- 
lite LRGs is not capable of reproducing the small-scale 
correlation function, no matter what concentration we 
use. We therefore relax the NFW assumption by allow- 
ing the inner slope of the profile to vary. Recall that the 
NFW profile has a logarithmic slope dlnp/dlnr of -1 at 
scales much less than the scale radius , and -3 at scales 
much larger than r^. We assume a new density profile 
for satellite LRG galaxies that is similar to NFW, except 
that the inner slope is no longer fixed to -1, but is a new 
free parameter —7: 

Ps 



p{r) 



3-7' 



(7) 



This reduces to NFW for 7 = 1. A model of this form 
has been used by papers that stud y the inner slope of the 
dark matter den sity profile (e.g.. iFukushige et al.ll200j : 
IReed et al.l[2005h . 

In order to use this new profile in our modeling, we 
need to compute the pair distributions Fcsif) and Fss{r), 
as described in § 13.31 While Fcs{r) is the profile itself, 
Fss(r) is the convolution of the profile with itself and 
quite non-trivial to calculate analytically for arbitrary 
values of 7 and concentration. We therefore calculate 
Fssir) in a numerical fashion. We make a dense grid 
of 7 and concentration values, and at each grid point 
we create an artificial spherical halo by putting down 
30k particles that satisfy the radial profile for that grid 
point. We then measure the pair distribution by count- 
ing all the particle pairs in our constructed halo. Once 
we have a table of Fss{r) functions on our grid, we can 
estimate Fgs (r) for any values of 7 and concentration by 
interpolating in the grid. 

As before, we wish to keep the same number of free pa- 
rameters in order to more fairly compare different mod- 
els. We thus keep Mq fixed to Mmin and we keep /gai 
fixed to unity. Therefore, we now vary the following 4 
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free parameters: Mmin, Mi, a, and 7, and we refer to this 
model as PNMG. We find a best-fit model with a reduced 
of 0.82 (x^ of 4.11 with 5 degrees of freedom). The 
dashed-dotted red curves in Figure [1] show this best fit. 
The PNMG model is clearly successful in fitting the M06 
small-scale data. Allowing the inner slope of the satellite 
LRG density profile to become steeper than is exactly 
what was needed to match the data. The value of 7 is 
well constrained and our MCMC yields 7 = 2.06 ± 0.21. 
It is certainly not surprising that the inner slope of the 
satellite density profile is similar to the slope of ^(r) at 
small scales because, as we argued in § 13.31 most LRG 
pairs should be central-satellite pairs whose pair distri- 
bution Fcs{r) is essentially the density profile itself. 

We have established that neither P{N\M), nor the con- 
centration of satellite LRGs, are sufficient to explain the 
M06 small-scale data, and that a profile other than NEW 
is needed. We thus now allow our model to have more 
than 4 free parameters and we vary both 7 and /gai- 
Since our density profile is no longer NEW, there is no 
reason to keep the concentration fixed to what was found 
for NEW dark matter halos. In our final model, we thus 
vary 5 parameters: Mmin, Mi, a, /gai, and 7, and we re- 
fer to this model as PNMCG. Our goal for investigating 
this model is to determine exactly what constraints the 
M06 data place on the density profile of satellite LRGs. 

We find a best-fit model with a reduced oi 0.71 
(x^ of 2.82 with 4 degrees of freedom). Figure O shows 
our l-,2- and 3-ct contours for 7 and /gai. As before, we 
find that the satellite LRG profile is much steeper than 
NEW (for NEW, 7 = 1 and /gai = 1), with 7 = -2.17 ± 
0.12. As we argued in ^ 13.11 most of the LRG satellites 
reside in halos of mass close to Mi. At our best-fit value 
for Ml (10^'*'^^/i~^Mo), the dark matter concentration 
fitting formula from § 13.11 gives a halo concentration of 
c = 4.25. Applying our 1-a range for /gai implies that 
Cgai ~ 0.1 — 4.1. Concentration values of unity or less 
mean that the scale radius is larger than the virial 
radius, which essentially means that the density profile 
retains its inner slope most of the way out. In other 
words, our fit shows that the density profile for LRG 
satellites is consistent with a simple isothermal profile. 

5. DISCUSSION 

Our results show that the distribution of satellite 
LRGs within dark matter halos requires a steeper in- 
ner density profile than NEW, which suggests that these 
galaxies are poor tracers of the dark matter distri- 
bution at these scales. The density profile of dark 
matter halos in the ACDM model has been measured 
extensively using high resolution N-body simulations, 
and recent inner profile measure ments seem to con- 
firm the NEW predictions ( Navarro et"al] I2008D 

- with so me slight deviations found in separate work 
fPicmand et all 120041: IFukushige et "all l2004t IReed et all 
I2b05t ,Del Popolo fc Kroupall2009D . However. NFW does 
not consider baryons, which can affect the dark mat- 
ter density profile at small scales. The interaction be- 
tween baryons and dark matter is addressed by the 
adiabatic contraction model that describes the gravita- 
tional effect of baryons on dark matter as the gas con- 
denses and sinks to the center of the dark matter po- 
tential well. The gravitational influence of the baryons 
draws the dark matter in, and this can steepen the 
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Fig. 2. — Allowed density profiles for LRG satellite galaxies 
within their dark matter halos. Shaded regions show the l-tr {red), 
2-CT (blue) and 3-cr (green) allowed regions for 7 and /g^i for the 
PNMCG model. Filled circles mark the locations of the NFW pro- 
file (7 = 1, /gal = 1) and the Singular Isothermal Sphere (SIS) 
profile (7 = 2, /gal ^ 0). 



density pr ofile (iGne din et al.' '2004^; ' Romano- Dfaz et aLl 
2008; Wein berg et a l. 2008; Sommcr-Larscn & Limousi^ 
Although it has been shown that the inn er pro - 
file can significantly steepen ( Gustafsson et all 12006( 1 . 
the majority of results show only a moderate steepen- 
ing. This has also been obs ervationally confirmed us - 
ing galaxy-galaxy lensing by iMandelbaum et aT] ()2006D 
who find that the mass density profile of LRG clusters 
is consistent with NFW. LRG satellites therefore have a 
steeper density profile than dark matter even with the 
effects of baryons taken into consideration. 

It is not necessarily surprising that LRGs are poor 
tracers of the dark matter density distribution within 
halos. LRGs presumably live in subhalos, which can cer- 
tainly have a different distribution than their host ha- 
los. .Nag ai fc Kravtsov. (2005) found that subhalos actu- 
ally have a shallower profile than dark matter at larger 
scales, but this has not been studied for the massive ha- 
los and small scales we consider here. It is difficult to 
model this regime because simulations must have, both 
a very large volume that contains many cluster-sized ha- 
los, and high mass and spatial resolution to resolve scales 
of ~ lO/i^^kpc. In fact, on these very small scales, it is 
likely that dark matter subhalos have already been dis- 
rupted by tidal forces, while the LRGs, being smaller and 
denser, have survived. So pure dark matter simulations 
may be insufficient to predict these LRG results. How- 
ever, semi-analytic models that include galaxies and can 
have arbitrarily high resolution should be able to make 
these predictions fsee lKitzbichler fc Whitell2008l for semi- 
analytic modeling of the very-small scale clustering of 
lower luminosity galaxies). 

In any case, our results have implications for the mod- 
eling of LRG clustering because a standard NEW profile 
cannot describe the spatial distribution of LRGs on scales 
< 0.03/i^^Mpc. It would be interesting to see if lower lu- 
minosity galaxies exhibit the same behavior or if this is 
simply a feature for LRGs. It would also be interesting 
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to see if LRGs maintain their steep density profile at high We thank Joanna Dunkley and David Hogg for useful 

redshift. discussions and comments. We thank Zheng Zheng for 

providing us with his best-fit 2-halo term for LRGs. 
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